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Abstract 

The algebraic matrix hierarchy approach based on affine Lie sl(n) algebras leads 
to a variety of 1 + 1 soliton equations. By varying the rank of the underlying sl(n) 
algebra as well as its gradation in the affine setting, one encompasses the set of the 
soliton equations of the constrained KP hierarchy. 

The soliton solutions are then obtained as elements of the orbits of the dressing 
transformations constructed in terms of representations of the vertex operators of the 
affine sl(n) algebras realized in the unconventional gradations. Such soliton solutions 
exhibit non-trivial dependence on the KdV (odd) time flows and KP (odd and even) 
time flows which distinguishes them from the conventional structure of the Darboux- 
Backlund Wronskian solutions of the constrained KP hierarchy. 



1 Introduction. The Algebraic cKP Model 



A large class of 1 + 1 soliton equations belongs to the so-called constrained KP (cKP) hier- 
archy. Some of the most prominent members of this group are the KdV and the nonlinear 
Schrodinger equations of the AKNS model. The cKP evolution equations possess the famil- 
iar Lax pair representations with generally pseudo-differential Lax operators which emerge 
naturally as reductions of the complete KP hierarchy Lax operators (PI)- Conventionally, 
the constrained KP hierarchy is obtained from the KP hierarchy by a process of reduction 
involving the so-called eigenfunctions. The eigenfunctions appear in the constraint relations 
introducing a functional dependence between initially infinitely many coefficients of the KP 
Lax operator. This scheme results in the pseudo-differential cKP Lax operator of the type 
C = L K+ i + Y^iLi ^id^^h where L K+ i is the differential operator of the K + 1-order, while 
are the eigenfunctions of C. In general, C possesses a finite number of coefficients 
which enter the soliton equations and depend on all (ti,t2,t 3 , . . .) isospectral time flows of 
the KP hierarchy. In the special M = case with the cKP Lax operator being a purely 
differential operator £ = Lk+i we encounter dependence on only some of the original time- 
flows of the KP hierarchy. The most simple example (K — 1, M — 0) is the KdV hierarchy 
with only odd time flows present. 

The soliton solutions for the cKP models have been found in (||) for the arbitrary K 
and M = 1 case using the Darboux-Backlund technique. Generalization to an arbitrary M 
is simple and was given in (||) (see also (Q)). These solutions appeared in the Wronskian 
form in terms of the eigenfunctions of the "undressed" C = d K+1 Lax operator. 

Here, we will present an alternative algebraic viewpoint of the constrained KP hierarchy. 
In this setting the algebraic dressing methods will provide new soliton solutions which appear 
to differ from the conventional form of the Darboux-Backlund Wronskian solutions due to 
a non-trivial mixing of the KdV-like versus KP time flows. This will be shown explicitly in 
the example characterized by K = M = 1. 

In an algebraic approach to the constrained KP hierarchy (H) the soliton evolution 
equations emerge as integrability conditions of the following matrix eigenvalue problem: 

d 

L^ = ; L = (D - A- E) ; D = I— (1) 

ox 

with the matrix Lax operator L belonging to Kac-Moody algebra Q = sl(M + K + 1). The 
integrable hierarchy is determined by the choice of gradation of Q. By varying the Kac- 
Moody algebras together with their gradations one is able to reproduce from the matrix 
hierarchy of eq. ([3]) the nonlinear evolution equations of the cKP hierarchy. 

We will be working with a simple setting in which the matrix E in eq.(|l]) has gradation 
1 with respect to gradation specified by the vector (@): 

s=(l,0 1 _^,L__ ; _ I l) (2) 

M K 

We call this gradation an intermediate gradation as it interpolates between the principal 
Sprincipai = (1, 1, • • • , 1) and the homogeneous one s homogeneous = (1, . . . , 0). As it is well- 
known the Wilson-Drinfeld-Sokolov (0 ||, ||, 0, |ll|]) procedure gives, respectively, the (m- 



)KdV (@) and AKNS (@||) hierarchies in these two limits. 



1 



Alternatively, the gradation s can be specified by an operator: 

K 

Q s = J2 X M+a-H° + {K+l)d (3) 
0=1 

Here Xj are the fundamental weights and d is the standard loop algebra derivation. Corre- 
spondingly, E stands for 

k ^ 

E = H E a M +a + E -{c lM+1 +-+ aM+K ) ( 4 ) 

0=1 

It is a non-regular (for M > 0) and semisimple element of Q. 

The matrix A in eq. (|l|) contains the dynamical variables of the model. A is such that it 
has gradation zero and is parametrized in terms of the dynamical variables % r«, U a and v 
as follows: 

M K 

A = Y, (<li p i + r i p -i) + E U M+a (a M+a ■ # (0) ) + v c (5) 

i=l a=l 

where P+; = Ev} n , , , >, , i = 1, 2, . . . , M, and c is a central element of Q. 



2 The Heisenberg Subalgebra and the Vertex Operator 

For a regular element E in the conventional Drinfeld-Sokolov approach the isospectral flows 
are associated with the Heisenberg subalgebra which can be identified with Ker (ad E). Here, 
due to non- regularity of E the Heisenberg algebra is associated to the center of Ker (ad E). 
It consists of the following three separate sets of operators: 

1) A homogeneous part of sl(M) , i — 1, 2, . . . , M — 1 

«y = ^ ,aW ; n,,^T) (6) 

2) A principal part of sl(K + 1) , a = l,2,...,K 

K+l-a a 

A a = -4- p( n+1 } (7) 

^ i o+n(Jf+l) - C/ a i+ Af+ai+M+l+--+«i+Af+a-i ~T~ -(a i+M +ai+M+i+-+ai+M+K-a) ^ ' 

8=1 8=1 

3) "A border term" 



aO iW+A tl , x JV / iW 



These relations provide a parametrization of the Heisenberg subalgebra in terms of ele- 
ments: 

bN,a = *4.Ar =a+n (^ + i) , &jV,0 = ^N=n(K+l) 5 ^N,i = (9) 
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where a = 1, 2, . . . , K, i = 1, 2, . . . , M — 1. The Heisenberg subalgebra elements from 
enter the oscillator algebra relations (we put c = 1): 



i>N,a 




= N5N+N'8a,K- 


_i_6 ; a, 6 = 1,2,.. 


.,K 


bN,0 


, bN',o} 








bN,i 


, bN',j ] 


= NSn+n 1 5ij 


; i,j = 1,2,..., M 


- 1 



(10) 

(11) 

(12) 



Define next the Fubini-Veneziano operators: 



OC 



V-iri) n oo AO -n(K+l) 

r^i Z „ , s v - ^n(K+l) z 



Qi<i<M-i(z) = — ; QM(z) = iJ2 " ^ , n — (1 3 ) 

n=l U n=l U \ JX L > 

oo J\a z a+n(K+l) 

QM +a (z) = , E a+ni " +1 ^^ u ; a=l,2,...,K (14) 

a + n(K + l) 

The corresponding conjugated Fubini-Veneziano operators Q*(z) are obtained from (|X3|)- ([14D 
by taking into consideration rules JC^\ = JC\~ n) , (A° n(K+1) ) = A°_ n(K+1) , (A a a+n{K+1) ) = 

In reference (flliJl) we found the step operators of sl(M + K + 1) associated with the 
Cartan subalgebra defined by the Heisenberg subalgebra (|5|). That in turn enabled us to 
find the corresponding simple root structure for sl(M+K+l) with the intermediate grading. 

Knowledge of roots and the Fubini-Veneziano operators is all what is needed to write 
down a compact expression for the general vertex operator in the normal ordered form: 

V a (z) = z^f 1 ''* 1 ' exp (id* ■ Q\z) ) exp (id ■ q) exp (a ■ pin z) exp (id ■ Q(z)^j (15) 



with M + f^-component root vector d described in ([{14)) and M + i^-component vector Q 
having components described in (|I~3|)-(|T4]). The zero-mode vectors p and q have only first M 
components different from zero according to (p)i = pi Q{M — i) and (q)i = qi 8(M — i). They 
satisfy relations [pi , qj ] — —iSy. Furthermore, Pm is equal to v4° =0 from expression @. 
An explicit example of the s/(3) vertex will be given below in section 4. 



3 The Dressing Technique and the Tau-function 

The dressing technique (||15||) deals with reproducing of the nontrivial part E + A of the Lax 
matrix operator from eq. ([!]) by the gauge transformations involving generators of positive 
and negative gradings applied to the semisimple element E: 

E + A = Q E Q' 1 + (d x Q) Q- 1 (16) 
E + A = (b^t) E (T^B) + (dxB^r) (t^B) (17) 



where B~ Y T contains positive terms and is an expansion in the terms of negative grading 
such that 6 = exp (j2i<o = 1 + ^ ^ + • • •• From expressions ( JIB]) and QT7D we obtain 
two alternative formulas for the same term A of grade 0: 

A = - \E, 6 { - 1] ] or A = -B- 1 (d x B) (18) 
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The term 0( ^ of grade —1 can be expanded as 

M+K 

fl(-i)- V e^E {0) + 

a=M+l 

M M 

+ E ^^(U..^) + E ^&U* + * (is) 
«=i 1=1 

where we included all possible terms of grade —1 according to 
Therefore 

M+K 

E,e^}= J2 t\^ (0) + C(-K + i + ... + % +K )-ff (0, + c) 

a=M+l 
M 

+ E^ e(«A/ + i,-«,-.-«M + i) -E_( Qi+ ... +QM) 

M 

+ E^ £(-aAf+X---«Ar+K>aj+-+aAf+K) ^q,+...+Qa/ (^0) 

Comparing the last expression with the field content of A as given by ([5]) we obtain relations 
for expansions parameters used in (|T9|): 

v = -9^ x) , U a = -0<T 1) + e£ 1) , r l = -6t 1) e(* M+1 ,-a l -...-aM +l ) 

Ql = —0\ e(-a A f+i-— -aiM+K ,ai+— +<xm+k) (21) 



We now work with representation of A as given in eq. flTq). We split the grade zero element 
B in a product B = B\B 2 with i?x containing the grade zero sl(M) elements and 



K 



B 2 = exp £ M+O a M+a • # (0) +P-C (22) 

0=1 

Accordingly, eq. (|T8|) becomes A = —B 2 ~ 1 B{ 1 ( y d x Bi)B2 — -B^c^f^) and A can be rewritten 

as 

K 

A=-J2 d x <f>M+a a M+a ■ # (0) -d xP -c + 0{sl{M)) (23) 

a=l 

where 0(sl(M)) contains all possible terms belonging to the sl(M) algebra. 
Comparison with (||) yields 

U M + a = -d x 4>M +a ; v = -d x p (24) 

We define a family of the first order differential matrix operators Cn = d/dt^ — An, 
N = 1, . . .. The hierarchy is then formulated in terms of the zero-curvature equations for 
the Lax operators 

[C N ,C M }=0 (25) 
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expressing commutativity of the higher time flows. The zero-curvature equations imply the 
pure gauge solutions for the potentials A N : 

C N = V^-V- 1 (26) 
ot N 



The starting point of the dressing method ([|15[]) is the vacuum solution v = U = r = q = 0. 
The corresponding Z/ vac ) = D — E matrix Lax operator together with higher flows operators 
Cn ,N > 1 for the vacuum solutions are expected to be recovered via (|26| ) from which 
is expressed entirely by the Heisenberg algebra associated with the center of Ker (ad E). 
Explicitly, for our model 

* = ^ (vac) = exp ( f] t N b {N) ) (27) 

\N=1 J 

with 6at given in (ffl). 

We define the tau-function vectors as: 

|r > = ¥ vac ^ h |A ) ; \r M+ a) = ^ (vac) h ¥ vac ^ \X M+a ) (28) 

They are associated with the constant group element h and the highest-weight vectors 
| Ao) , IAm+o) such that 

a M+a ■ # (0) |A > = ; a M+b ■ H<®\\ M+a ) = 5 a , b \X M + a ) ; a, b = 1, . . . , K. (29) 

Assuming that h allows the "Gauss" decomposition of \|>( vac ) h \I>( vac ) 1 in positive, nega- 
tive and zero grade elements we get for the tau-function vectors from (p8|) an alternative 
expression: 

|r ) = e-^lAo) ; \r M+a ) = Q- 1 B- 1 \X M+a ) 

0-1 = ^(vac) h ^(vac)- 1 ^ . g-1 = ^(vac) h ^(vac)- 1 ^ ^ 

for a — 1, . . . , K. Like before, we make a splitting B^ 1 = B^B^ 1 in ( |30|) and notice that 
B^IXm+o) = I^Af+a) for -Bi being an exponential of sl(M) generators. Inserting B 2 from 
(2^) we find 

|r ) =0- 1 |Ao)e ( ^ ) ! |r M +a) = © _1 S" 1 |A M+a > e ("^+-) (31) 



r (0) = exp (-p) ; rjJL = exp (-p - (f) M+a ) (32) 



Denote 

1 M+a 

Accordingly, expanding 9 _1 as below eq.(|H]) we find 

\tm+o) 



r (o) 

T M+a 



l-9^-...)\X M+a ) (33) 



and similarly for \tq)/tq°\ We find by comparing with relation ([21]) that 

T\ = e(a M+1 ,-a l -...-a M+1 ) ( A^f + i | -E^+...+aA/+i I / T M+1 (^4) 

Ql = e(-a M +l---a M +K ,a l +...+a M+K ) {X Q \E^ ai+ _ +aM+K ) |to)/t q ° (35) 

U M+a = -d x In (rf/rSU ; 1/ = -0, In (r (0) ) (36) 
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The multi-soliton tau functions are defined in terms of the constant group elements h 
which are the product of exponentials of eigenvectors of the Heisenberg subalgebra elements 

h = e Fl e F * ■■■e Fn , [b N , F k ] = J$ F k , fc = 1, 2, . . . , n . (37) 

As seen from eq. (|37j ) for such group elements the dependence of the tau-vectors upon the 
times t N can be made quite explicit 

n (k) 

\r a ) = J] exp(e£>" ** Fk ) | Aq ) (38) 
fc=i 

The multi-soliton solutions are conveniently obtained in terms of representations of the 
"vertex operator" type where the corresponding eigenvectors are nilpotent. 



4 The sl(3) Example: Solitons of the Yaijma-Oikawa 
Hierarchy 



We apply the above method to the particular case of sl(3) with M — K — 1. From eq.( 
the surviving elements of the Heisenberg subalgebra are in this case: 



h (2n+l) 



>N=2n+l,a- 



(n+1) 



-a 2 



b (2n) = b N=2n>0 = VsXi-H^-^=6 n , Q 



b (k) b (k') 



(39) 
(40) 



kSk+k' for both even and 



and they satisfy the usual Heisenberg subalgebra 
odd k. 

The structure of eigenvectors of Heisenberg subalgebra facilitates construction of mul- 
tisoliton solutions according to fl57p and (|3"%|). In the current example we find that the 
eigenvectors and their corresponding eigenvalues (in notation of fl3T|)) are 



E, 



-2n- 



1 (Eg ~ E^) + 

_ 2z 2n+l 

E &1+&2 = ^11^? + 



(2n+l) _ z 2n+l 



■ i , (2n) 



E 



H (n) _ £ A 



n. 



(2n+l) 
U a2 



-2n-l 1 . , ,(2n+l) _ . 



2n+l . (2n) 



(42) 
\/3z 2r (43) 



We now realize the above eigenvectors by the nilpotent vertex operators. The construction 
involves the Fubini-Veneziano operators defined in terms of the Heisenberg elements as in 
eqs. (HD-©: 



b {2n+l) z - 



2n-l 



«6Z 



2n + 1 



(52(2) = ? - ip\nz + z ^ 



2n~ 



■2/1 



(44) 
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where the zero mode momentum p = = ■ — c/2a/3 satisfies [q , p] = i. The 

step operators from (|4l| ) - (|43|) are then realized, from the algebra point of view, as vertex 
operators via: 

#fii <"> E {i,V3)(z) = ^ z3/2 ■ ex P (iQi(*) + iV3Q 2 (z)) : (45) 



E &2 » E^ t0) (z) = ~:exp{-2iQi{z))--e i ** (46) 

E &1+&2 <-> E^^z) = V2z 3 ? 2 : exp (-iQ^z) + iy/3Q 2 (z)) : (47) 

and similarly for the negative root step operators, with change of sign of % in exponentials. 
A care has to be exercised in applying this correspondence within the setting of the Fock 
space with the vacuum vector being | A2) , since (Ao|-Ea 2 l^o) = — 1/2 while (A 2 |-Ea 2 |^2) = 1/2 
as seen from expression fl4"2|). Similar consideration applies for the products of E a /s vertex 
operators like £V x _^E,_ 1 ^ which produce E^ 2 ,o)- 
Introduce the notation: 

V CiA (z) = z d ' /2 : exp(ic i Q 1 (z)+id i Q 2 (z)) : (48) 

It is not difficult to establish the following correlation function: 



Xr. 



n*r 2VTi 2 n u 1 -^ [(, i -, i )^+^)]^ /2 (49) 

3=1 l<i<j<n \ Zl Z 3 J 

for cr = 0, 2 and with 

00 00 

r C j,dj (^j) = C i ^2n+l^| n+1 + X/ 4/ ^2n^| n (50) 
n=0 n=l 

When substituting V^.^^) by E Cj)dj (zj) one encounters extra phases originating from the 
Klein factor in eq.((46|) and from the character of the [ A2) vacuum as discussed below eq. 
(|47|). The latter gives rise to the factor exp (im/2) Y%=i c j) f°r the A 2 correlation function 
as verified on several examples. 

Recall that for the problem in hand the Lax matrix operator from ([!]) with A from (||) 
and E from @) specifies to 

^ 

1 - vc (51) 
-U 2 J 

where A is the usual loop parameter. In terms of the tau- vectors we have from (|34|) -(|36"1) the 
following n-soliton representation of the components of the Lax operator 

1 1 n 

r = ^<A 2 |i^ + Jr 2 > = ^ <A 2 ^ (52) 
r 2 r 2 j=i 

1 1 n 

9 = ^(Ao|^- Q Jr ) = ^(Aoi^l_ a2 ^ ) n(l + ^^(%))^ r IV (53) 

U 2 = -^ln(r (0) /r 2 (0) ) ; v = -d x In (r (0) ) (54) 




where 



-(o) 



n 



(55) 



3=1 



Using association between the step operators (fH|)-(|43|) and the vertex operators (f4"5|)-(f4?|) 
we can rewrite the step operators appearing in (|5^) and (|53"D as 

1 



E, 



(0) 



2lTX 



(56) 



We now calculate the zero-curvature equations (|25|) [D — E — A , d tn — A n ] = for the 
first two non-trivial cases of n = 2, 3. We expand A n = J27=o A n {i) where the index i in the 
parenthesis equals grading with respect to Q s = X 2 ■ H° + 2d. We choose A 3 (3) = E$ +E^l 2 
and A 2 (2) = v 7 ^ • in order to ensure truncation of the expansion. This method yields 
for the first non-trivial case (n = 2) the evolution equations 



= d t r + d 2 r + rd x U 2 — qr 2 — U 2 r ; = d t U 2 + d x (rq) 
= d t q- d 2 x q + qd x U 2 + rq 2 + U 2 q 

where we defined for simplicity t = \f?>t 2 . The evolution equations for t 3 are 

t 








4 



dx {rd x q - qd x r) 



3 
2' 



3 



d t ,r - dlr - ^d x rd x U 2 - -jd 2 x U 2 + ^rU 2 d x U 2 + ^U 2 d x r 



3 
— i 
2 



3 2t T O 3 2 r\ 

+ -qr U 2 + -rqd x r - -r d x q 



d t3 q - dlq + \d x qd x U 2 + \qd 2 x U 2 + \q\J 2 d x U 2 + hj 2 2 d x q 
-g r(7 2 + -^rqd x q - -q d x r 



(57) 
(58) 



(59) 



(60) 



(61) 



(£)V n/2) , £ 



These equations follow also from the conventional Sato equations d tn C - 
applied to the scalar cKP Lax operator C = (d — U 2 )(d + U 2 — qd~ l r). We note here 
that the simple reduction of the matrix Lax operator from ([51]) yields the scalar spectral 
problem C\\ = A% with the scalar Lax operator L\ — (d + U 2 )(d — U 2 — rd~ 1 q). Both 
scalar Lax operators are related by a conjugation and the Darboux-Backlund transformation: 

A = (d + u 2 )c*(d + u 2 y\ 

We now present few examples of the soliton solutions (p2])-(^5]) satisfying the above 
evolution equations. 

1) One-soliton solution, n = 1. With h = (1 + E^ 2j o(zi)) we recover the standard m-KdV 
one-soliton configuration with r = q = and 



-(0) 



,(-2x Zl -2t 3 zi 3 ) 



r W = 1+ - e (- a .*-aW) (62) 

2) Two-soliton solutions, n = 2. For h — (1 + ^-2,0(21)) (1 + -^i ^(^2)) we find Tq ' 1 and r 2 ' ) 
as in (R3) but now with g 7^ and equal to 



q 



V2 



22 e 



(*3 z 2 3 +i Z2 2 +a; z 2 ) 



1 + 1 (-2.ft-.2W) f^±^)W 



«1 - z 2 



(63) 
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while r = 0. Similarly, for h = (1 + £/_ 2 ,o(^i))(l + -E^ -^3(22)) we find r^O but q — 0. 
For /i = (1 + ^i >v ^(«i))(l + ^l.-VaO^)) we find tlmt botn 9 7^ and r ^ 0: 

_ l+cr/2 2-cr/2 ( I Zi+t«l 2 +t3«f+X Z 2 -t22 2 +t3Z2 ) 

r (0) = 1 + ( _ 1) W 2 ) 2 ^ gg e _ L ; . = 0,2 (64) 

(zi-z 2 )(zi + z 2 ) 2 

■s/2z 2 e*- ~* 2 2 2 + :c - z 2+t322 ) y/2 z% e( tZl2+XZl+t3Z i ) 



r (0) ' y r (o) 

r 2 r 



(65) 



3) Three-soliton solutions, n = 3. As an example we take here h = (1 + £'_2,o(^i))(l + 
^i > ^2))(l + S 1| _V3(^)). We find 



-(0) 



= 1 _|_ (_1)W2) i e (-2x2 1 -2t 3 2i 3 ) 
~ 2-0-/2 _ l+cr/2 

_|_ ( 1]W 2 ) 2 e (t3Z2 3 +tZ 2 2 +XZ 2 +t 3 Z3 i -tZ 3 2 +XZ 3 ) /gg\ 

(z 2 -23) (^2 + 2:3) 2 

_|_ 2: 3 2 ~ CT/2 ^2 1+g/2 + Z 2 ) + Z 3 ) ^ _ 2 x Zi _ 2 13 Zi s +t Z2 2 +x Z2+f3 Z2 3_ t zs 2 +x Z3+t3 23 3 ) 

( z 2 ~ Z 3 ) ( Z 2 + Z 3 ) 2 ( Z 1 - Z 2 )( Z 1 - Z 3 ) 



V2 Z3 ei-^ +xM ( x + 1 /r (o) (6?) 

(68) 



- 23 



In the above examples U 2 and f can be obtained from fl54|). We note that we only kept the 
explicit time dependence on times t n with n < 3, which was enough to verify the evolution 
equations fl5"7|)-(|6"l~D. 

The novel feature of the above soliton solutions is that they mix exponentials exp fe^=i tnZj 
which represent a typical time dependance for the KP solutions with pure KdV-like time 
dependance of the type exp (j2^ =0 t 2n +i^ 2n+1 ) involving only odd times. 

The exception is provided by the pure KP type of solution in eq.(^)- (|65[) , which can 
also be obtained as a Wronskian arising from the Darboux -Backlund transformations. 

5 Conclusions 

We presented here new soliton solutions for the cKP hierarchy. They exhibit a nontrivial 
mixing of KP times t n with all indices n's and KdV-like times t 2n -i with only odd indices. 
From the point of view of intermediate gradation, such a mixture is indeed very natural. In 
fact, we are able to obtain pure KdV solutions, pure KP solutions (meaning all times) and 
the arbitrary mixtures thereof, just by varying the constant group elements h of the dressing 
orbit. 

We were not able to obtain these solutions from the conventional Darboux-Backlund 
method involving Wronskian representation. Recall, that such Wronskian is given in terms 
of the eigenfunctions fi satisfying d n fi = d x f\ for all n. This type of time dependence fits 
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only a subclass of our soliton solutions, namely those which are of the pure KP type (meaning 
that all t n 's with all n's are present without the nontrivial mixture with KdV times). Based 
on studying many examples of these solutions, we believe that the Wronskian method works 
equally well in these cases. 
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